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ABSTRACT 


TMs  report  discuses  *  definition  of  the  effective  coupling  coefficient  for  the 
body-force  class  of  electromechanical  transducers.  This  definition  is  based  on 
the  distribution  of  energy  in  the  transducer  and  is  therefore  particularly  applica¬ 
ble  to  the  'distributed  parameter'  type  of  transducer.  Furthermore,  the  definition 
results  in  a  frequency  dependent  coupling  coefficient  which  is  valid  for  all  modes 
of  vibration. 

The  definition  is  applied  to  some  common  types  of  longitudinal  vibrator  and 
the  results  are  compared  with  ex  peri  menu  for  one  particular  device  (a  partially 
excited  stack  or  ceramic  discs).  The  Influence  of  parasitic  elements  on  the  system 
coupling  is  examined. 
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DYNAMIC  COUPLING  COEFFICIENTS 
FOR  DISTRIBUTED  PARAMETER  PIEZOELECTRIC  TRANSDUCERS 


INTRODUC  TIOK 

Distributed  parameter  transducers  may  be  described  with  some 
success  by  analogous  circuits  containing  frequency-dependent  circuit 
elements,  but  the  extraction  of  an  effective  coupling  coefficient  from 
such  a  circuit  is  usually  done  with  approximation  techniques  wb'di  have 
doubtful  validity  for  devices  wth  high  material  coupling. 

This  report  exploits  a  definition  of  the  effective  coupling  based  on 
ti  e  distribution  of  energy  in  the  vibrating  system  and  applies  it  to  vari¬ 
ous  specific  transducers.  This  definition  not  only  produces  exact  values 
of  the  coupling  factor  for  all  frequencies  and  modes;  it  also  sheds  light 
on  the  somewhat  enigmatic  "meaning"  of  the  coupling  coefficient  con¬ 
cept  for  composite  transducers. 

In  this  report,  the  development  begins  with  a  derivation  of  the  equa¬ 
tions  of  motion,  followed  by  a  complete  electromechanical  description 
of  three  types  of  distributed  parameter  transducers  (all  of  the  "longi¬ 
tudinal  vibrator"  class).  This  work,  of  course,  is  available  in  many 
places  in  the  literature;  it  is  reproduced  here  to  establish  nomenclature 
for  the  main  body  of  the  report. 

The  third  section  introduces  the  definition  of  coupling  to  be  used 
and  applies  it  to  the  three  transducers  described  in  the  preceding  sec¬ 
tion.  A  general  discussion  of  coui  ling  coefficients  for  composite  trans¬ 
ducers  follows.  A  theoretical  description  of  a  partially  excited  longi¬ 
tudinal  vibrator  is  compared  with  experimental  data  in  the  final  section. 

The  transducers  considered  in  the  report  are  all  piezoelectrically 
excited,  but  the  concept  applies  equally  well  to  the  magnetostrictivc 
devic  ?s. 


DERIVATION. OF  THE  LONGITUDINAL 
VIBRATOR  EQUATIONS 


THE  CONSTANT  FIELD 
LONGITUDINAL  VIBRATOR 

Consider  a  segmented  longitudinal  vibrator  of  length  2l  which 
vibrates  in  the  x-direction  symmetrically  about  its  center.  The 
lateral  dimensions  of  the  bar  are  small  compared  to  t  so  that  lateral 
vibrations  are  ignored  and  the  shape  of  the  r ross- sectional  area.  A. 
is  arbitrary  (rectangular,  circular,  and  annular  shapes  are  common). 
The  deflection  £(x)  of  the  plane  at  position  x  ,  and  the  stress  at 
this  point,  T(x),  obey  the  differential  equation 

,•*£(*)  ,  HQ  -  0  (1) 

ds 


for  harmonic  excitation  at  angular  frequency  u 
material  is  p  . 


The  density  of  the 


Numerous  electrodes  placed  in  the  y -/  plane  and  evenly  spaced 
with  separation  t  along  the  x-axis  produce  an  electric  field  whose 
x  component  is  S(x>  and  an  electric  display  ":i.ent  with  x  compo¬ 
nent  D(x).  The  electrodes  are  connected  in  parallel  in  the  usual 
fashion.  These  electric  variables  are  connected  with  th.-  stress  and  the 
strain,  S(x),  by  the  following  piezoelectric  equations  of  state: 

JS(.>  -  (2) 

(D<»)  -  d,,T<«)  +  (*)  - 


The  strain  is  the  gradient  of  the  deflection,  S-d£/dx.  Assuming  that 
the  electrodes  are  spaced  so  closely  that  the  electric  f:eld  strength 
does  not  vary  along  the  length  of  the  material  (d8/d*«0),  Eq.  (1)  be¬ 
comes  the  familiar  Helmholtz  Equation: 

-  a.  (3) 

ox 


The  general  solution  of  Eq.  (3)  is 


£(*)  «  Bj  sis  k*  +  Bj  co*  lex  ,  (4) 

where  the  wave  number  k«t»/c*  .  and  the  speed  of  sound  in  the  mate¬ 
rial  at  constant  field  is  cj  -(psjjir*.  Applying  the  boundary  conditions 
on  f  : 

f<”  '  °  )  <5> 

f  («)  -  j 

results  in  the  exact  solution  of  the  differential  equation, 

m  -  r,  •  161 

*  itakl 

and,  immediately,  the  distribution  of  strain: 


S(«) 


d(  t,  com  kx 
dx  “  C-  mm  k( 


(7) 
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Use  of  the  first  equation  of  state  now  gives  the  stress  distribution  in 
the  vibratoz: 


T(s) 


cos  k» 

si*  kt 


4  8  1 

*»  •  I 


(8) 


(The  subscript  on  £  is  a  reminder  that  the  field  is  now  regarded  as 
a  spatial  constant. ) 


If  the  ends  of  the  sample  are  free  lT(t)-0],  Eq.  (8)  establishes 
the  relation  between  the  two  amplitude  parameters,  and  Sa:- 

k  col  kt  -  .  (9) 


The  resonant  frequencies  of  the  longitudinal  vibrator  depend  on  ihe 
electrical  terminations.  If  the  terminals  of  the  transducer  are  short- 
circuited,  the  natural  frequencies  are  found  from  Eq.  (9)  with  ,=  0  : 


iwiMi&MHaia 


co«  k*  i  »  0 

*B*  -  -5; 


•  3i  5( 


(The  superscript  E  indicates  constant  voltage  conditions.)  The  cur¬ 
rent  flow  between  two  electrodes  is 

it  -  j<uADa  (11) 


where  Da  is  u.e  value  of  the  displacement  current  in  the  region  be¬ 
tween  x.  and  xa_,  .  The  total  current  flow  is  the  sum  of  the  ia's: 


I  -  j«AlD  -  A 


•’"*(')  jjf 


(0  °- 


D(s)<k. 


Under  the  many-electrode  assumption  (t/*»i),  D(s)  may  be  regarded 
as  a  continuous  function  of  x.  It  is  obtained  by  inserting  Eq.  (8) 
in  the  second  equation  of  state: 


»  fii  £  k  .  €t  A  d«  \p 

^  +  »  r77)  •' 

**  \  •»»**»/ 


Recognizing  the  material  coupling  coefficient 

k. 

^  a  •  _  


”  «*  eT 

•irii 


and  performing  the  integration  of  Eq.  (13)  results  in  the  total  current 


1  -  >"7  j^r  ^.+€» 


Under  open-circuit  conditions,  there  is  no  current  flow.  Putting  Eq.  (16) 
equal  to  zero  yields  the  open-circuit  deflection/field  relationship  which 
must  coexist  with  Eq.  (9).  Combining  these  two  produces  the  equation 
for  the  allowed  frequencies  at  constant  current: 

k 


k*  I  co*  k‘{  - - 


l  -  k2 
*  *»» 
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The  open>circuit.  voltage  is  simply 


V 


8  t 


(18) 


Examination  of  Eqs.  (10)  and  (17)  shows  that  the  short-circuited 
frequencies  (constant  voltage  resonances)  depend  only  on  elastic  con¬ 
stants  and  dimensions,  while  the  open-circuited  frequencies  (constant 
voltage  anti-resonances)  depend  also  on  the  piezoelectric  properties  of 
the  material  as  indicated  by  the  presence  of  kM  in  Eq.  (17).  Further 
more,  if  one  measures  the  fundamental  (m  =  1)  resonance  and  anti¬ 
resonance  of  an  unloaded  stack,  the  material  coupling  may  be  found  by 
forming  the  ratio  r  •  ai1/*)*  and  noting  that 


k*l  -  fk*»  >  rl. 

2 


(19) 


Now  Eq.  (17)  becomes* 


«■  r  -  cot  r  - 

1-kJ,  2  2 


(20) 


The  electrical  admittance  of  the  device  is  found  by  dividing  the 
current  in  Eq.  (16)  by  the  voltage,  8#t: 


Y.  -  ~  -  >-A 

c  * 


[ *»  kt 


(2») 


Because  this  analysis  considers  a  lossless  system,  the  electrical  admit 
tance  is  entirely  imaginary.  Rewriting  Y  in  the  form 


f  A6”1 


i  kt 


kl 


A(l-kJ,> 


’) 


l??'. 


shows  that  the  admittance  has  the  form  of  two  suscepiances  in  parallel. 
The  second  term  represents  the  susceptance  of  a  capacitor  of  value 


At  €* 

C.  * 


(23) 


'"This  equation  is  •  "ell-known  re  lit  ion  between  material  coupling  and  the  resonam-to-anti- 
resooant  frequency  spread.  See:  BerUncoun  et  al. ,  froc.  Institute  of  Radio  Engineers,  vol.  48,  p. 


[ 
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while  the  first  term  is  called' the  motional  susceptance.  By  introducing 
the  characteristic  mechanical  impedance  of  the  material,  Zc  =  pcA, 
this  motional  susceptance  takes  the  form 

B  _  (Ad>1/«*>t)>  WJ  (24) 

-jZ#«tkl  "  -JZt  ootkl  4 

where  M>Adta/s*,t  is  called  electromechanical  transformation  ratio. 
An  equivalent  circuit  which  has  the  admittance  Y>  at  its  terminals 
is  shown  in  Pig.  1. 


Fig.  X  -  CpwiM  Circuit  of  the  Constant  Field  Lan^iudtiul  VI beat oc 

Examination  of  Eq.  (21)  shows  that  the  resonant  frequencies  de¬ 
fined  bv  Eq.  (10)  —  the  constant  voltage  resonances  —  cause  Y#  -  •  , 
while  those  frequencies  defined  by  Eq.  (17)  —  the  constant  current 
resonances  —  produce  Y#  •  0. 

The  capacitance  C.  is  called  the  blocked  capacitance.  It  may  be 
interpreted  as  the  capacity  of  the  segmented  bar  when  it  is  clamped  to 
prevent  motion;  th-  dielectric  constant  applicable  for  this  calculation 
is  designated  "longitudinally  clamped":  (I  -kaj). 

The  analysis  of  a  thin  bar  in  longitudinal  vibration  due  to  transverse 
electric  fields  ("31"  mode  of  operation)  follows  the  form  of  the  pre¬ 
ceding  theory  closely;  the  "31"  subscript  replaces  the  "33"  on  all 
material  constants.  The  reason  for  the  similarity  is  that  both  methods 
of  excitation  share  the  common  characteristic  that  the  electrically  pro¬ 
duced  strain  is  uniform  in  the  direction  of  motion.  This  is  achieved  in 
t£e  "33"  case  by  using  closely  spaced  electrodes  placed  normal  to  the 
e  vector  and  i  ;  the  "31"  case  by  using  continuous  electrodes  along 
the  sides  of  the  bar. 


If,  however,  a  bai  vibrates  oo  that  an  appreciable  portion  cf  the 
acoustic  wavelength  fits  between  the  electrodes,  there  will  be  a  variation 


of  field  in  this  interval.  This  may  occur  for  higher  modes  of  the 
"33"  vibrator  just  discussed  (large  m),  or  for  the  fundamental  mode 
of  a  column  of  ceramic  with  electrodes  only  on  the  ends  (*  -  +!)  and 
none  within  the  material  (i.  e. ,  a  solid,  rather  than  segmented,  stack). 
Such  electromechanical  behavior  is  characterized  by  d8/dx  4  0 
(spatially  variable  field)  and  is  described  in  the  next  part  of  the  report. 

THE  VARIABLE  FIELD 
LONGITUDINAL  VIBRATOR 

The  analysis  of  this  type  of  vibration  parallels  that  just  discussed. 
The  same  geometry  and  coordinate  system  are  used,  but  it  is  appro¬ 
priate  to  employ  a  different  set  of  Equations  of  State: 

S(x)  -  .J.TWtlyDCi) 

(25) 

£<*>  -  -«„T(x>  +  0j,L(:  . 

It  is  interesting  to  see  how  the  material  constants  with  superscript  D 
are  related  to  those  with  superscript  E.  For  constant  D  (or  D  =  0), 
Eq.  (2)  becomes 

•  <26> 


Eliminating  2  from  this  pair  yields 


This  is  Hooke' s  Law  for  the  constant  D  case;  the  elastic  constant  is 
smaller  than  s'j  and  is 

kj,)-  (28) 

The  alectric  constant  pT  appearing  in  Eq.  (25)  is  the  reciprocal 
of  the  permittivity  of  the  material,  eT  .  The  piezoelectric  constant  g 
is  related  to  d  by 


4. 


The  first  of  Eq.  (25)  is  inserted  in  the  general  differential  equation, 
Eq.  (1).  The  resulting  equation  is  simplified  by  recognizing  that, 
although  8  may  have  spatial  variations  along  the  x  axis  of  the  bar, 
the  electric  displacement  must  be  constant  and  thus  dV/dx  -  0 . 

The  two  electrical  conditions  are,  as  before,  constant  voltage  and 
constant  current.  The  first  relation  is  obtained  by  making  the  total 
voltage  developed  ac-oss  the  two  electrodes  vanish: 

E  m  J*  8  (z)  dx  -  0  .  (30) 

The  open-circuit  conditions  are  established  by  making  D  =  0.  The 
two  sets  of  natural  frequencies  that  result  are  found  from: 


short  circuited:  <t*  I  cot(*Et}  » 


open  circuited: 


cot  (**  i  )  «  0 


The  wave  numbers  **  and  are  different  from  those  in  the  preceding 
section  since  different  elastic  constants  are  used  for  this  type  of  vibra¬ 
tion: 

(«*  -  o*/Cj  ,  ft*  -  W*/Cj  ). 

Note  that  the  constant  voltage  resonances  depend  on  coupling  and  the 
constant  current  frequencies  do  not.  This  is  the  reverse  of  the  behavior 
of  the  segmented  bar  described  in  the  preceding  section. 

A  circuit  which  has  the  same  electrical  impedance  as  the  vibrator 
is  shown  in  Fig.  2. 


-  j  Z0  cot/ci 


Fig.  2  -  Equivalent  Clrcui:  of  the  Variable  Field  LiHigitudinal  Vibrator 


The  values  of  the  circuit  elements  are 


Ae 


C  - 


11 
2 1 


(1-kJ,). 


N  . 


a  a. 


2.^1 


The  frequencies  which  satisfy  Eqs.  (31)  and  (32)  make  the  admittance 
of  this  circuit  infinite  and  zero,  respectively.  Equation  (31)  shows 
that  the  constant  voltage  resonances  are  displaced  below  their  values 
in  the  dB/o*  »  0  case;  the  overtones  are  not  integral  multiples  of  the 
fundamental. 


This  is  represented  graphically  in  Fig.  3  where  (tamct)  and  (a?  /kjs) 
are  plotted  against  *1 ;  the  intersections  of  the  two  curves  represent 
the  natural  frequencies. 


Note  that  the  roots  are  displaced  from  the  lower  modes,2  but 
approach  even  overtones  for  large  (*l). 

Furthermore,  for  low  values  cf  the  material  coupling  coefficient 
k3J  ,  the  slope  of  the  straight  line  increases  rapidly,  and  the  f  requeue  v 
disolacements  become  small  even  for  the  lower  modes. 

JThu  behavior  occur*  also  foi  shear  modes.  See;.  H.  F.  Tlersien,  journal  of  ihe  Acousncai 
Society  of  Amcnca,  vol.  36,  p.  53.  January  iSo3. 


4,  •-  c-v  ^ 


By  measuring  the  u*  and  iJ  for  the  fundamental  mode  and  forming 
the  ratio  r  =  u//u>*  ,  one  can  determine  the  material  coupling  from 

‘1,  ■  =r  “■  (t)  •  <33' 

This  equation,  and  its  companion  for  the  segmented  bar  Eq.  (20), 
assume  a  more  familiar  form  by  using  the  notation  f*  =  f*  +  Af.  Then3 

j  Seated  B*  TTr-if^T?)1 “(57?)  .  <34> 

I  (dB/dx-0)  1-kJ3  2  \  f  /  \  f  / 


En  6-elec  tioded  Bex. 

OD/dx  -0) 


■”  •  i  (‘-7) “(i  7) 


These  relations  are  illustrated  in  Fig.  4. 


a  (CO.  20) 
•-  f£0-  33) 
C.  f*G.  63) 


rig.  4  -  Ceramic  Material  Coupling  Coefficients  as  Deiermmed  from  the  Ratio  of  Aiui-Resonam 
to  Resonant  Frequencies  for.  A.  Segmented  Sat  and  S<de-Electroded  Bar  toe  6z-0), 

B.  Enri-Elecuoded  Bai  (40  -  0>  ;  and  C.  Effective  Coupling  for  Lumped  Constant  Transducer 


Berltncourt,  op.cit. (footnote  1). 


.0 


Equation(34)  may  also  be  used  to  find  kj,  for  samples  with  their 
electrodes  parallel  to  the  displacement.  This  is  because  d&/dx  -0 
is  satisfied  for  side-electroded  bars  also. 


THE  PARTIALLY  EXCITED 
LONGITUDINAL  VIBRATOR 

The  type  of  vibration  discussed  in  this  section  results  if  the  seg¬ 
mented  bar  is  driven  with  a  given  number  of  the  electrodes  on  either 
end  disconnected  so  that  only  the  central  portion  of  the  bar  (symmetri¬ 
cally  located  about  the  origin)  is  electrically  excited. 

Clearly  the  ceramic's  elastic  properties  are  different  in  the  two 
regions;  therefore,  the  wave  equation  must  be  solved  separately  in 
each  portion  of  the  bar  and  the  additional  constants  eliminated  by 
joining  the  two  solutions  at  the  point  of  discontinuity. 

Let  the  parameter  y  designate  the  fraction  of  the  bar  that  is 
excited  electrically  (0<y.<l),  and  let  the  subscripts  "d”  and  "f" 
designate  the  driven  portion  (|s|  .<  yi  )  and  the  free  portion 
(y!  <|*|  <  !',  respectively. The  equations  of  state,  Eq.  (2),  apply  to 
die  d-section  and  those  of  Eq.  (25)  to  the  f-section.  In  using  Eq.  (25) 
it  should  be  recalled  that 

D(m)  -  0  far  t  >  s  >yt  , 

since  those  electrodes  are  disconnected  from  each  other.  If  the  elec¬ 
trodes  in  the  f- region  were  tied  together  but  not  connected  to  the 
electrical  generators  different  sort  of  motion  would  result,  but  that 
case  is  not  considered  here. 

The  analysis  in  the  d-section  follows  the  derivation  of  the  constant 
field  longitudinal  vibrator  up  to  the  equation  for  the  displacement  func¬ 
tion,  Eq.  (4).  Making  the  vibration  symmetrical  about  the  origin  with 
fd(0)  -  0  and  specifying  the  displacement  at  the  plane  of  discontinuity, 
£d(y* )  -  4y  ,  yields  the  particular  solution 


where  the  wave  number  K4  applies  only  to  the  d- region. 


and 


*4*  * 


(37) 


(38) 


Differentiating  Eq.  (36)  produces  the  strain  function;  this  is  then 
combined  with  the  first  equation  of  state  to  get  the  stress  distribution 
in  the  d- section  of  the  bar: 


T4(») 


1 


cat  K4s 


(39) 


A  procedure  similar  to  the  above  is  used  to  find  6f(s)  and  Tf(s) . 
The  boundary  conditions  used  are: 

a.  Matching  displacements  at  the  y  interface:  £f(yl)  -  £y  - 

b.  Unloading  the  free  end:  Tf(t )  -  0 . 


The  results  are 


ff(*> 


co<  A(  co*  K,i  ♦  «io 
cot  Af  co«  y  Af  +  sin  y  Af 


and 


Tf(s) 


Kf£y  coi  K(i  -  co:  A(  ii*  K(< 
cos  y  Af(cot  Af  ♦  taa  y  Af) 


(4.) 


(41) 


where  the  wave  number  is  w/c^  Af  *-  •  The  remaining 

condition  on  the  two  solutions  is  that  of  matching  the  stresses  at  s  -  yl  . 
Since  there  are  no  more  undetermined  constants  in  the  problem,  this 
operation  will  produce  a  fixed  relation  between  the  two  amplitude- 
determining  parameters,  sad  The  resulting  equation  is  the 

analog  of  Eq.  (9)  for  the  partially  excited  bar: 


X 


/  C,  \  ^ 

-v-rA  I  ■;.  • 


The  natural  frequencies  of  vibration  when  the  electric  terminals 
are  short-circuited  are  those  which  satisfy  Eq.  (42)  when  »  0 . 
Designating  the  allowed  values  of  A4  for  this  condition  by  the  super  - 
script  E  (constant  voltage),  and  using  the  relation  Af  -  A4  Vl  -k|, 
and  the  trigonometric  identity 


cot  y  A  -  cot  A 
cot  y  A  cot  A  +  1 


.(i-y)  A, 


one  obtains  the  following  implicit  relation  for  A*  as  a  function  of  y 


(tsafl-yWl-k*,  A*',  MoyAj-Vt^kj,  . 


The  constant  current  resonances  are  found  by  making  Dd*  »  0. 

The  distribution  of  electric  displacement  along  the  bar  is  found  from 
the  second  equation  of  state: 

D4<»>  *  +  ) 


K  .  cos  K  .z 

.T  D  i.  il%  «  v  4  < 

-  £,»t.(l-kM)  + - w -  - — T  ’ 

*■  •  * 3  -*  *u»  y  A . 

*n  •  I 


Performing  the  integration  and  equating  the  result  to  zero  yields  th~ 
open-circuit  voltage/deflecllou  i  elation: 


f*J  C«  +  ty  “ 


This  result  is  combined  with  the  fixed  Cy/S*  relation  given  in 
Eq.  (42)  to  obtain  the  open-circuit  resonances. 


’  Ai  | 


cot  y  - 


VT-fci 

’  11 


>  U -  y)  Vl-kf 


IS 


(46) 


Equations  (43)  and  (46)  reduce  to  Eqs.  (10)  and  (17),  respectively, 
when  y  -  1.  Combining  Eqs.  (43)  and  (46)  to  obtain  an  expression 
for  kM  in  terms  of  the  resonant  to  anti-resonant  frequency  spread  is 
a  more  formidable  chore  than  it  was  for  the  fully  excited  cases  treated 
earlier.  Furthermore,  it  is  doubtful  that  ih“  result  of  such  work  would 
be  very  useful  since  it  is  unlikely  that  knowledge  of  ku  would  be 
sought  with  experiments  on  partially  excited  bars,  it  is  instructive, 
however,  to  observe  the  variations  of  w1  and  «“  with  y  .  The 
values  of  these  two  frequencies  for  y  =  1  serve  as  convenient 
normalizing  values  and  one  may  examine  the  dimensionless  ratios 


and 


«“(y) 


(47) 


(48) 


With  this  notation  and  Eqs.  (10)  and  (17),  the  y  -dependent  A 
variables  become 


Aj(y)  -  (»/2)  8*  (y) 


(49) 


A^(y)  -  r(»/2)«‘(y) 


(50) 


The  solutions  of  Eqs.  (43)  and  (46)  for  5 *(y)  and  tff y)  appear  in 
Fig.  5  for  the  specific  cas'  k„  =  0.64  (r  =  1.22).  The  constant 
current  solution  is  displaced  upward  by  the  factor  r  on  the  graph 
so  that  the  curves  refer  all  relative  f regencies  to  -*(!). 

Toulis  has  published  data  from  an  experiment  involving  a  partially 
excited  stack4  with  an  r  =  1.  22  for  y  =  1.  These  data  are 


V  J.  Toulii,  "Redefinition  of  Effective  Bectiomechamctl -Coupling  Factor,"  Journal  of  ihe 
Acouaucal  Society  of  Amenca,  35,  no.  12,  p.  2024,  November  1963. 
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superimposed  on  die  curves  of  Fig.  5.  The  failure  of  the  experimental 
anti-resonant  frequencies  to  coincide  with  the  theoretical  curve  may  be 
due  to  the  fact  that  the  measurement  of  anti- resonant  frequencies  ;s 
generally  more  difficult  than  resonant  frequencies  and  mechanical 
losses  (which  were  ignored  in  the  analysis)  are  more  influential  at  ws. 


THE  ELECTROMECHANICAL  COUPLING  COEFFICIENT 
THE  QUASI-STATIC  CASE 

In  the  introductory  sections  of  Hersh's  excellent  report  on  coupling5 * 7 
there  are  listed  for  review  manyfamiliar  definitions  of  the  coupling  coeffi¬ 
cient.  Many  of  these  definitions  concentrate  on  extending  the  coupling 
concept  from  electrical  circuit  analysis  to  the  realm  of  electromechani¬ 
cal  devices.  Such  techniques  are  indeed  useful  since  the  description  of 
electromechanical  devices  through  the  use  of  equivalent  circuits  is  a 
popular  method  of  analysis.  In  the  opinion  of  this  writer,  however,  these 
extensions  lose  much  of  their  usefulness  when  applied  to  distributed 
parameter  transducers,  since  the  circuits  which  describe  such  devices 
often  contain  frequency-dependent  elements. 

This  objection  causes  one  to  examine  another  class  of  definitions:, 
those  which  seek  to  describe  coupling  in  terms  of  the  distribution  of 
electrical  and  mechanical  energies  in  the  transducer.  One  definition 
from  this  group  which  has  particular  applicability  to  the  distributed 
parameter  transducers  examined  in  this  report  is  that  described  by 
Vigoureux  and  Booth*  Skudrzyk,*  and  Beckmann*  This  view  of  coupling 
is  also  exploited  by  Berlincourt  et  ah  in  their  contribution  to  the  recent 
book.  Physical  Acoustics.* 

The  distinguishing  feature  of  this  approach  is  that  it  considers  the 
total  energy  of  the  transducer  system  to  be  the  sum  of  three  terms: 
the  purely  electrical  and  purely  mechanical  energies  without  coupling, 
and  the  mutual  (or  shared)  energy  resulting  from  coupling.  The  square 
of  the  coupling  coefficient  for  a  particular  mode  is  proportional  to  the 
square  of  the  mutual  energy  divided  by  the  product  of  the  two  self  energies. 

That  this  delinition  follows  easily  from  the  linear  equations  of  state 
for  piezoelectric  transducers  is  shown  by  Beckmann  in  the  following  manner- 


5  John  F.  Hersh,  Coupling  Ci  efficients.  Acoustics  Research  Laboratory,  Harvard  University, 
Technical  Memorandum  no.  40,  pp.  4-25,  1957. 

*  P.  Vigoureux  and  C.  F.  Booth,  Quarts  Vibrators  and  Their  ApplicaUops.  H.  M.  Stationery 
Office,  London,  p.  58,  1950. 

7  E.  Skudrzyk,  Pie  Grundlagea  det  Akuttlk,  Springer-Verla^  Wien,  p.  485,  1954. 

*  R.  Beckmann,  'Some  Applications  of  the  Linear  Piezoelectric  Equations  of  State,  ~ 

I.R.E.  Transaction:  oc  Ultrasonics  Engineering,  PGUE-3,  pp.  66-59,  May  1955. 

*  W.  P.  Mason,  Ed. ,  Physical  Acourtics,  vol.  »,  Pan  A,  Academic  Press.  New  York, 
p.  189  ff.,  1964. 
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Select  a  pair  of  equations  of  state  and  write  them  in  the  form 

(W5(,V£i.»e- 


(51) 


The  mechanical  energy  density  in  the  volume  where  these  equations 
apply  is 


liSvT? -iSs?  T*Tu  +  i2d 
2  A  A  2  A  1*  2  • 


iA  'A  ■ 


+  ®»a* 


(52) 


while  the  electrical  energy  density  is 


ilD.e*  -  is  8* T„+ij 


p  p 


Pp  p 


P*  P  _ 
i  '“P  m 


(53) 


where  the  asterisk  indicates  a  complex  conjugate,  W2  represent 
the  purely  elastic  energy,  W2  the  purely  electric  energy,  and  i\12 
the  mutual  piezoelectric  energy.  The  square  of  the  electromechanical 
coupling  coefficient  is 


k* 


(54) 


Equations  (52)  and  (53)  are  quite  complicated  in  the  general  form; 
however,  if  mechanical  and  electrical  boundary  conditions  for  a 
particular  type  of  excitation  leave  only  a  few  of  the  ST  and  CD  products 
non-zero,  the  coupling  coefficient  may  assume  a  more  simple  form. 


This  may  be  illustrated  by  considering  a  sample  with  the  applied 
field  parallel  to  the  axis  of  polarization  (6^,1  0,  Sj-8,-0)  and  all 
stresses  normal  to  this  axis  zero  (T,  =  T2  -0,  T,  ^  0).  Under 
quasi-static  conditions,  T2  and  Sj  are  uniform  throughout  the 
piece  andthe  applicable  equations  of  state  Eq.  (51)  reduce  to 


t-q.  U): 


Tj  +  eJl  ''a 


(55) 


The  energy  densities  of  Eqs.  (52)  and  (53)  are  now  quite  simple,  and  one 
finds  the  total  to  be 


( 


1 

+  - 
2 


i }  j 


<“>6) 


I  ■  i  i*  T!  +  2 

2  *33  1 


The  above  procedure  applies  to  calculations  using  the  so-called 
homogeneous  equations  of  state  in  which  the  dependent  and  inde¬ 
pendent  state  variables  are  each  intensive  or  extensive,  respectively. 
As  Beckmann  points  out,*  if  one  uses  the  "mixed"  equations  cf 
state,  the  ratio  of  energies  Fq.  (54)  is  not  k  but  rather  is 
(k*/(l  -  k1 ).  This  important  point  is  summarized  in  Table  1. 

Title  1 

APPROPRIATE  RELATIONS  BETWEEN  COUPLING 
COEFFICIENT  AND  ENERGY  RATIO  FOR 
VARIOUS  EQUATIONS  OF  STATE 
(The  wlniin  t  hAcsti  the  tnaapoic  of  • 


10  Beckmann,  op.  clu,  p.  *6. 
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DYNAMIC  COUPLING  COEFFICIENTS 


For  non-static  conditions  the  energy  densities  Wt,  W  j .  W  u  are 
integrated  over  the  volume  of  the  transducer  to  account  for  the  distribu¬ 
tions  of  stress  and  field  which  occur  in  the  distributed  parameter  cases. 
The  total  energies  thus  obtained  are  then  used  to  form  the  ratio  in 
Eq.  (54). 

This  procedure  will  be  applied  to  the  three  examples  oi  such  trans¬ 
ducers  which  have  been  considered  in  this  report. 

In  the  case  of  the  constant  field  ’ongitudinal  vibrator  (the  seg¬ 
mented  bar)  one  has  for  the  strain  energy,  with  piezoelectric  effects 
ignored. 


»,  “5«?,/T,dv- 


For  T  one  may  use  Eq.  (8)  with  d  „  =0.  Calling 
the  stress  at  the  center  of  the  bar,  Eq.  (58)  becomes 


l  T2kJ 

[  »  -»*  A  ■■■  ■■  ■  I  cot2  kx  dx 

1  2  **  ui’k t  J 


1  *»AkT.* 


[let  +  'in  let  cot  kt  ] 


Since  the  field  is  constant, the  electric  energy  is  simply 


-  -  :T,  /V  av  « I  At  e2 . 
2  2  *  2  13 
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The  mutual  energy  is 


"5a«/T8dv 
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Now  Eq.  (54)  yields  the  "dynamic"  (or  "effective")  coupling  coefficient 

m  *«  _  **»» _ 2  sis*  kl _ 

*  "  t  m  "  eT  ,»  kf  (kl  +  si*kl  cos  kl ) 


11  |  kl  (kl  +  sis  kt 


cos  kt) 


This  effective  coupling  is  plotted  against  the  mode  number,  m  =  kt/(ir/2) 
for  kj,  =0.64  in  Fig.  6  (solid  curve).  Note  that  the  coupling  goes 


kj  -  kj3  F(k!) 
k?  k* 


- -1- - iL.  FOcI) 
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2  sin2  kl 

kl  (kl  +  sin  kl  coskl) 


k.3  -  0.M 


1.0  2.0  3.0 

■n-kt  /{»/2) 

Fig.  6  ■  Frequency -Dependent  Effective  Coupling  Coefficient*  foe  Two  Type*  of 
Longitudinal  Visitors 
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to  zero  for  even  m  since  these  are  prohibited  modes  of  vibration. 11 
At  zero  frequency  the  effective  coupling  is  equal  to  the  material 
coupling,  ku,  since  the  stress  distribution  becomes  constant. 

Of  particular  interest  is  k,  at  the  fundamental  resonance: 


Since  there  is  a  constant  p rope -‘tonality  between  these  two  factors 
one  could  obtain  k.  as  a  function  of  the  frequency  spread  ratio  r 
by  multiplying  Curve  A  in  Fig.  4  by  (•/•*)*  X  0.90.  When  this  is  done, 
however,  the  resulting  curve  differs  only  slightly  from  the  function 
k  -  y/T^T7tr ,  This  latter  expression  is  the  effective  coupling  of  a 
simple  lumped  constant  equivalent  circuit  for  a  one-degree-of -freedom 
transducer;  it  is  usually  seen  in  the  form 

kj-l-<f,/f/.  (63) 


The  k#  of  Eq.  (63)  is  plotted  on  Fig.  4  as  Curve  C. 


For  the  constant  D  case  treated  previously  (the  end-electrodedbar), 
a  similar  procedure  is  used  to  find  the  dynamic  coupling  factor.  The 
principal  difference  is  that  the  mixed  equations  of  state,  Eq.  (25),  are 
used  in  the  analysis,  and  therefore  Eq.  (54)  does  not  apply.  The  energy 
densities  in  an  mlinitesimal  volume  element  are  lound  by  multiplying 
the  equations  of  state  by  T  and  D,  respectively. 


w,  +  *sa 


_  V  +  W 
12  T  2 


(64) 


Following  the  rules  in 


Table  1 
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(65) 


"This  c»n  alio  be  seen  in  Equation  (9):  For  k!  «  m(r/2)  and  ■  even,  the  volt  age 
requited  to  maintain  vibration  becomes  infinite. 


i 


And  for  the  dynamic  case,  of  course. 


I  k l  ~  D*  dV] (/ sjj  T*  dV] 

m  kl>  [JDTavl* 

'  1-kJ,  l /D*  dV)  [  /T*  dV) 


(66) 


The  computation  to  evaluate  this  expression  for  the  end-electroded  bar 
turns  out  to  be  exactly  the  same  form  as  that  preceding  Eq.  (62),  and 
the  result  is 


kl  kj,  j  > 

1-kJ  1  -  kjj  I  (x  I  +  sia  *  t  cos  *  t )  J  ((>7) 


The  results  appear  graphically  as  the  dashed  curve  in  Fig.  6  for 
k„  =  0.  64.  Note  that  as  w  0  and  the  stress  distribution  in  the 
bar  flattens  out,  the  expression  in  brackets  has  a  limit  of  unity  and 
again  k.  (static)  =  k„ .  Since  the  resonant  frequency  of  the 
dD/dx-0  bar  is  coupling-dependent  (see  Eq.  (31)),  it  is  difficult 
to  evaluate  Eq.  (67)  at  resonance.  The  anti-resonant  frequencies, 
though,  are  coupling -invariant  and  thus 


1  -k. 


1-k 


* t 
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Unlike  the  previous  case  there  is  no  constant  relationship  for  k(/k,$ 


For  a  specific  k„  ,  however,  one  can  evaluate  k#  versus  r 
by  applying  the  above  equation  to  Curve  B  in  Fig.  4.  The  result  is 
again  graphically  indistinguishable  from  Curve  C.  One  may  conclude, 
therefore,  that  Eq.  (63)  is  valid  for  the  experimental  determination 
of  the  effective  coupling  of  the  unloaded  end-electroded  bar  at  anti¬ 
resonance. 

This  close  approximation  t*  the  actual  effective  coupling  by  a  Simple 
formula  is  fortuitous.  The  reason  for  the  agreement  lies  simply  in  the 
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fact  that  the  effective  stiffness  of  the  transducer,  though  frequency- 
dependent,  changes  only  slightly  in  the  interval  from  resonance  to 
anti- resonance.  The  variable  r  in  Eq.  (63)  is  obtained  by  meas¬ 
urements  at  two  frequencies,  while  the  value  of  coupling  which  results 
is  associated  with  just  one  frequency.  In  lumped-constant  circuits  no 
difficulty  is  encountered  by  this  since  one  constant  value  of  coupling 
coefficient  applies  to  all  frequencies  for  which  the  circuit  10  valid  — 
the  distribution  of  state  variables  throughout  the  transducer  is  constant 
with  frequency.  When  such  constancy  is  not  true,  the  measurement  of 
the  ratio  of  two  separate  frequencies  may  not  be  an  adequate  means 
of  determining  the  coupling.  The  derivation  of  Eq.  (63)  rests  on  the 
assumption  that  the  mechanical  reactance  of  the  transducer  varies  in 
a  specific  manner  (i.  e. ,  that  of  a  series  L-C  circuit),  and  the  inclu¬ 
sion  of  reactances  which  vary  according  to  some  other  scheme  (e.g., 
a  transcendental  equation  such  as  that  in  Fig.  2 )  destroys  the  validity 
of  the  derivation. 

Although  this  discussion  shows  that  the  use  of  Eq.  (63)  will  not 
produce  exactly  the  correct  value  of  effective  coupling  for  some  com¬ 
posite  transducers,  the  equation  still  has  great  utility  for  comparative 
studies  of  a  specific  transducer.  That  is,  although  the  numerical  value 
of  k.  obtained  may  be  incorrect,  one  may  quickly  check  the  relative 
effect  of  design  modifications  on  the  coupling  coefficient  by  using 
Eq.  (63). 

One  further  point  of  interest  is  that  the  effective  coupling  of  the 
constant  and  variable  field  longitudinal  vibrators  diminishes  at  each 
ou«-^co»l»cly  higher  uveriuiit;  resonance,  inis  is  seen  clearly  in 
Eq.  (62).  for  increasing  (kl  )  the  numerator  maintains  the  periodicit 
while  the  denominator  causes  the  maxima  of  k.  to  approach  zero 
asymptotically.  This  phenomenon  is  due  to  the  numerous  nodes  in  the 
stress  p  'em  at  higher  modes.  It  is  evident  from  the  definition  given 
as  Eq.  (o4)  that  the  stress-field  product  should  not  vary  greatly 
throughout  the  body  in  order  to  maintain  high  coupling. 

This  effect  can  sometimes  be  used  to  increase  the  effective  coupling 
Regions  with  full  electrical  excitation  but  low  stress  levels  may  be 
removed  from  the  transducer.  For  instance,  in  flexural  bars  and  plates 
the  low  stresses  near  the  neutral  plane  deteriorate  the  coupling.14  If 

lIR.  S.  Woollen,  "Theory  of  the  Piezoelectric  Flexural  Diik  Transducer  with  Applications 
to  Underwater  Sound,"  USL  Research  Report  no.  490,  S-F001  0i  04-1,  pp.  71-74,  1960, 


these  layers  of  ceramic  are  replaced  by  low-modulus  inert  material^ 
the  coupling  may  increase. 

DYNAMIC  COUPLING  OF  THE 
PARTIALLY  EXCITED  BAR 

According  to  the  definition  of  coupling  used  in  this  report,  vibrations 
which  do  not  directly  participate  in  the  electromechanical  conversion 
process  (i.  e. ,  are  not  electrically  active)  and  static  electrical  elements 
which  store  electrical  energy  without  being  coupled  to  the  vibrations 
through  the  piezoelectric  effect  (i.  e. ,  external  capacitances)  both  tend 
to  decrease  the  effective  coupling  coefficient  of  a  transducer. 

One  example  of  this  effect  is  the  composite  transducer  in  which  an 
electromechanical  driver  excites  other  electrically  inert,  resonant,  or 
non-resonant  structures.  A  transducer  coupled  to  a  fluid-filled  cavity 
is  such  a  device;  the  effective  coupling  is  lowered  even  if  the  cavity  is 
tuned  to  have  rero  reactance  at  the  natural  resonance  of  the  driver. 

An  experiment  to  test  this  principle  has  been  performed  by  Toulis.** 
It  consists  of  measuring  the  resonant  and  anti-resonant  frequencies  of 
a  segmented  multi-electroded  stack  of  ceramic  for  varying  degrees  of 
electrical  excitation.  This  experiment  was  described  previously, and 
the  frequency  data  appears  in  Pig.  5. 

What  is  the  theoretical  dynamic  coupling  of  such  a  partially  excited 
structure  f •  the  energy-balance  viewpoint?  One  finds  immediately 
the  electrical  and  mutual  energies  (in  half  the  bar)  to  be 

Vr)-;«I.A f  K*  (68) 
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since  S(  =  0.  The  elastic  energy,  though,  has  two  distinct  parts:  a 
contribution  from  the  "d* 1  section  and  one  from  the  "f"  section: 


(70) 


By  Eq.  (54)  the  dynamic  coupling  of  the  composite  bar  is 


k*(y) 


+  w^u  +  u./iy 


(71) 


Comparing  Eq.  (71)  with  the  work  on  pages  19  and  20  reveals  that 
(W^ /W1U4)V‘  is  the  effective  coupling  of  a  fully  excited  constant 
field  bar  of  half  length  yt  ,  while  (1  +  U,  /Ud)  is  a  -orrection  term 
tv  account  for  the  electrically  inactive  end  sections  of  L  partially 
excited  bar. 

Computation  of  the  ratio  Uf  /Ud  at  system  resonance  as  a  function 
of  y  is  tedious.  The  general  form  (that  is,  at  all  frequencies)  is 


lij  *i»J  y  '4lu to*1  ♦(*»«!  \(cos  A' -stay  \(co*y  A()(l  -cix1  Af)-2co«  \f(ua2  A(-  sin2  >  Af)| 

_  g  (y) - — - - - - 

*  \  1-k  Jj  (»»ny  Af*  c»h  AjCo*)-  A()  (y  A^i  »wy  A4co»y  Atf) 

(72) 


To  find  the  effective  coupling  ac  the  fundamental  resonance  as  a  function 
of  y,  one  colves  Eq.  (43)  for  the  resonance  parameters  Aj(r), 
forms  the  functions 


A?  *  A; 


(73) 


and  evaluates  S(>)  from  Eq.  (72). 


The  relative  maximum  of  k#(y)  near  y  =  0.7  results  because 
the  mass  loading  of  the  electrically  active  central  part  by  the  unexcited 
end  sections  causes  the  stress  distribution  to  become  less  severe  and 
raises  the  coupling.  For  small  y  ,  however,  the  large  amounts  of 
strain  energy  in  the  end  sections  lowers  the  system  coupling  again. 

The  simple  formula  Eq.  ( 6S )  becomes,  in  the  notation  of  this 
problem, 


1  -(S*  t#)1 . 


Once  again  this  approximate  expression  is  sufficiently  accurate  for 
good  engineering  over  all  of  the  range  of  k3J  presently  available 
(less  than  1  percent  error  for  all  Y  and  all  kJ3  <  0.65). 

Applying  Eq.  (75)  to  Toulis'  data  (Fig.  5)  produces  the  experi¬ 
mental  coupling  versus  y  points  plotted  in  Fig.  8.  The  theoretical 
curve  for  kM  =0.  64  is  also  drawn  for  comparison  (solid  curve).  An 
earlier  study  of  a  partially  excited,  side-electroded  bar  was  made  by- 
Beckmann  and  Parsons.14  Although  their  treatment  was  simplified  by 
neglecting  the  change  of  elastic  properties  of  the  material  between  the 
active  and  inactive  regions,  the  general  shape  is  similar  (Fig.  8, 
dashed  line). 

•*0|— 
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Fig.  8  -  Dynamic  Coupling  Coefficient  of  the  Partially  Excited  Bar  vrrsrn  y 


1<R.  Beckmann  and  ?.  L.  Parioo,  Piezoelectric  it),  :i.  M.  Stationery  Office,  London, 
p.  296,  1957. 
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